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ABSTRACT 

Electron scattering induces a polarization in the cosmic microwave background (CMB) 
signal measured in the direction of a galaxy cluster due to the presence of a quadrupole 
component in the CMB temperature distribution. Measuring the polarization towards 
distant clusters provides the unique opportunity to observe the evolution of the CMB 
quadrupole at moderate redshifts, z ~ 0.5-3. We demonstrate that for the local clus- 
ter population the polarization degree will depend on the cluster celestial position. 
There are two extended regions in the sky, which are opposite to each other, where 
the polarization is maximal, ~ 0.1(t/0.02) /iK in the Rayleigh-Jeans part of the CMB 
spectrum (r being the Thomson optical depth across the cluster) exceeding the con- 
tribution from the cluster transverse peculiar motion if v t < 1300 km/s. One can hope 
to detect this small signal by measuring a large number of clusters, thereby effectively 
removing the systematic contribution from other polarization components produced 
in clusters. These polarization effects, which are of the order of (vt/c) 2 r, (i>t/c)r 2 and 
(kT e /m e c 2 )T 2 , as well as the polarization due to the CMB quadrupole, were previously 
calculated by Sunyaev and Zel'dovich for the Rayleigh-Jeans region. We fully confirm 
their earlier results and present exact frequency dependencies for all these effects. The 
polarization is considerably higher in the Wien region of the CMB spectrum. 



Key words: cosmic microwave background - cosmology: theory 
general - polarization. 



galaxies: clusters: 



1 INTRODUCTION 

The intensity of the cosmic microwave background radia- 
tion is distorted in the direction of galaxy clusters. There 



are two basic effects of this kind: the thermal effect (gun 



& Zel'd 



aev fc Zel'dovich 1972) and the kinematic effect (Sunyaev 



h 198C) The pnlariTatinn nf trip CMB should hf 



changed dllfi to the presence of '-.v flust e rs in the line 



of sight as we ll ( Sunyaev fc Zel'dovich 1980 ; Zel'dovich & 
Sunyaev 198C; [Sunyaev fc Zel'dovich 1981 ). In fact there 
are several separate polarization effects but since the ex- 
pected amplitude of each of them is very small, until re- 
cently there was little interest to them. However, we have 
observed such a significant progress in CMB measurements 
over the period passed since 1980, that the possibilty to 
detect a CMB-induced polarization signal from galaxy clus- 
ters seems no longer unrealistic. Furthermore, the CMB ex- 
periments planned fo r the next decade, which include the 



borne experiments, e.g. POLAR (Keating et al. 1998), will 
already reach a level of sensitivity which is not far from what 
is needed for such a detection. This drastic change in the sit- 
uation is reflected in the fact that several new papers (Audit 
& Simmons 1998; Itoh, Nozawa & Kohyama 1998; Hansen & 
Lilje 1999) discussing the polarization effects in clus ters have 
appeared over t he past year or so. Some authors (Audit & 
Simmons 1998; Itoh et al. 1998) have claimed the ne cessit y 
to revise the original results of Sunyaev & Zel'dovich (1980). 
This gave us a first motivation to return to this problem. 

Another reason for us to begin this study was the dis- 
covery by COBE (see, e.g. Kogut et al. 1996) of a quadrupole 
componen t in the CMB angular anisotropy. Zel'dovich & 
Sunyaev (1980) showed that the presence of a quadrupole 



space missions MAP ( http:/ /map.gsfc.nasa.gov), SPO rt on 



the International Space Station (Cortiglioni et al. 1998 ), and 
Planck Surveyor (tittp://astro. estec.esa.nl/Planck; De Zotti 
et al. 1999), and a number of ground-based and baloon- 



component should cause a change in the CMB polarization 
in the direction of a galaxy cluster. Already our first estimate 
based on the COBE measurement lead us to the rather sur- 
prising conclusion that this polarization effect should gen- 
erally be stronger than any of the other polarization effects 
relating to galaxy clusters. This fact makes the quadrupole- 
induced polarization potentially very important for cosmol- 
ogy. Indeed, measuring the microwave polarization signal 
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towards distant galaxy clusters opens the unique opportu- 
nity to effectively probe the CMB quadrupole anisotropies 
at intermediate redshifts, z ~ 0.5-3, i.e. from positions in 
the Universe different from ours, thus gi ving us a way to 
reduce the cosmic- variance uncertainty ( Kamionkowski & 
Loeb 1! 97). Furthermore, by observing the effect in clusters 



(or any large-scale ionized gas clouds) at large z we could be 
able to follow the evolution of the CMB quadrupole, which 
should be different for universes with different cosmological 
parameters. 

The structure of the paper is as follows. In Section 2 
we give the basic formulae which allow one to find the an- 
gular distribution of the polarized radiation resulted from 
scattering of a given initial angular distribution of unpolar- 
ized radiation by a single resting electron. In Section 3 we 
calculate the amplitude and distribution over the sky of the 
polarization imposed by the CMB quadrupole on the signal 
from clusters. We study other polarization effects connected 
with galaxy clusters in Section 4, which is followed by a Dis- 
cussion of possible ways to detect the cluster polarization 
due to the CMB quadrupole. 



2 POLARIZATION DUE TO SCATTERING OF 
ANISOTROPIC RADIATION BY A SINGLE 
ELECTRON 

In this section, which lays the basis for our further analysis, 
we basica lly re peat the original treatment of Zel'dovich & 
Sunyaev (1980), which is also de scribe d in detail in a review 
paper of Sunyaev & Zel'dovich (1981). 

We consider the situation that a single electron being at 
rest scatters initially unpolarized radiation. The frequency of 
a photon does not change due to scattering (one can safely 
neglect the small change due to recoil when dealing with 
the C MB), and we therefore can, for the moment, ignore 



the spectral dependence of the radiation and consider only 



its angular distribution. This distribution can, in general, 
be presented as the sum of components proportional to the 
Legendre polynomials P n (/i): 1, fl, (fi 2 — 1/3), etc., where fx 
is the cosine of the angles between the wave vector of the 
incident photon and one of a set of directions that uniquely 
define our particular angular distribution. The incident ra- 
diation is described by a usual set of Stokes parameters, (7, 
Q, U, V), where Q = U = V = 0. Our aim is to determine 
the corresponding quantities (to be marked with primes) for 
the scattered radiation. The fourth Stokes parameter, V, 
describes circular polarization which is not produced in a 
scattering and will, therefore, not be mentioned further. We 
define the Stokes parameters of the scattered radiation with 
respect to the plane that contains a) the specific direction 
about which a given radiation anisotropy component has an 
axial symmetry (see above), and b) the wave vector of the 
scattered photon, so Q > 0, U = corresponds to a lin- 
ear polarization in the direction perpendicular to this plane. 
With this choice of the reference frame, it is easy to show 
( Chandrasekhar 1950; Sunyaev & Zel'dovich 1981) that the 
scattered radiation will be 



f'M = 



i + £<V-D( 



I(n) dfj,, 



(1) 



QV) = / ^(1- M 2 ) (V - I) m dfi, U'(p') = 0.(2) 



One can see that the angular scattering function for the Q 
Stokes parameter is of strictly quadrupole form (with respect 
to /i). This fact is of importance for our further considera- 
tion, because it means that, due to the orthogonality con- 
dition for Legendre polynomials, only the quadrupole com- 
ponent, neither the dipole nor the higher harmonics, in the 
incident intensity distribution will contribute to Q' . 

Now let the isotropic and multipole components of the 
intensity distribution be Io, -Zi/i, ia(/* — 1/3), etc. Then 
upon integration in (§)-(§) we derive 



/V) 



Io + -h 



(3) 



Q'((i') = -Hl-V' 2 ), [/'(//) =0. (4) 

As expected, the isotropic component remains unchanged 
upon scattering. If one now considers a situation that the 
optical depth of the scattering medium is very small, t«1, 
then only single scatterings of photons will play role, and it 
is straightforward to modify formulae (|^), (Q), which were 
obtained for an individual electron, for this case: 



Io + To h { fI ' 2 -l 



= Am') 



10 



oo 

Ji^' + o.9/ 2 (y 2 - 1) + J « P "(m') 



-/ 2 (1- M ' 2 ), 0V) = 0. 



(5) 



(6) 



Equations (^), (^|) show that electron scattering tends to 
smooth anisotropies in the intensity distribution and to pro- 
duce a polarization on the order o f Q'/I' — 0.1r/2/Jo(l — 
fx' 2 ) (Zel'dovich & Sunyaev 1980; sunyaev & Zel'dovich 



1981) 



3 POLARIZATION DUE TO THE CMB 
QUADRUPOLE 

The COBE satellite provided the first measurement of the 
quadrupole component in the CMB angular distribution. Al- 
though the accuracy of this measurement is very low - its 
amplitude, Q rm3 , is said to lie in the range [4, 28] /iK at 95% 
confidence (Kogut ct al. 1996) - it is very likely that the de- 



tection of the quadrupole is real, because the presence of 
a similar quadrupole signal, Q r ms- P s ~ 15/xK, also follows 
from the extrapolation of the CM B power spectrum f rom 
higher multipole moments to I = 2 (Bennett et al. 1996). 

The CMB quadrupole should cause the microwave ra- 
diation from clusters of galaxies to be polarized. Below 
we present a formalism that allows one to calculate the 
quadrupole-induced polarization to be measured from a 
cluster with given celestial coordinates for a given CMB 
quadrupole. Rather as an example (because of the large un- 
certainty present) we find the distribution of the polarization 
signal over the sky as implied by the COBE measurement, 
ft will be straightforward to update our current, rough esti- 
mates as new, more accurate measurements (the first of this 
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kind is likely to be provided by MAP) become available. Al- 
though, the formulae below are directly applicable only for 
clusters which belong to the local population, i.e. for mea- 
suring the quadrupole that is observed at the current epoch 
from Earth, one can easily extend our approach on the de- 
termination of quadrupoles as seen from highly redshifted 
clusters. 

The COBE quadrupole is defined by 5 parameters, 
which are essentially coefficie nts at spherical second-order 
harmonics flKogut ct al. 1996[ ): Qi = (19.0 ± 7.4 ± 8.2) (jlK, 
Q 2 = (2.1 ± 2.5 ± 2.7) JjK, Q 3 = (8.9 ± 2.0 ± 2.5) /xK, 
Qi = (-10.4T8.0T4.3) fjK, and Q 5 = (11.7±7.3±10.4) fiK. 
Here, the first quoted errors are statistical, which are fol- 
lowed by the COBE team's estimates of the present sys- 
tematic errors. One can construct from the Qi parameters 
a tensor Q a/ s with the following components 



Qxx — Qa 2^ 1: Q yy — 

Qxy — Qyx — Q§i Qxz — Q zx 



(7) 



The quadrupole signal, Iq, can be calculated for any direc- 
tion as Iq = ^2,Q a 0T a ri3, where the Cartesian coordinates 
r a are expressed through the Galactic coordinates I and b, 
i.e. X = cos/ cos b, Y — sin I cos b, and Z = sin fa. The Q a p 
tensor possesses the property Q xx + Q yy + Q zz = 0. 

We can next find the eigenvalues Ai, A2, A3 (let Ai be 
the minimum and A2 the maximum among the three Ai) and 
the corresponding eigenvectors el, el, el for the Q a /3 tensor 
in a usual way: Q a pei = Ajej. On introducing spherical co- 
ordinates, 9 and <j) (with the polar angle 9 measured from 
el), we can, using the property Ai + A2 + A3 = 0, cast the 
quadrupole signal in a simple form: 



Iq = Ai sin 2 6 cos 2 4> + A2 sin 2 9 sin 2 (j> + A 3 cos 2 9 

= \ x ^sin 2 9 cos 2 4> — — ^ + A H ^sin 2 9 sin 2 <fi — i \ ,(8) 

where X x = An, — A3 and X y — A2 — A3. The RMS amplitude 
of the quadrupole signal is given by 



Qrms — Jg^s _ Aa-Ay + X y ) . 



(9) 



For the COBE specific values (as we use these only for il- 
lustration purposes, we will quote no errors for the numbers 
to follow), we find Ai = —25 /iK (minimum), A2 = 23 [IK 
(maximum), A 3 = 2 /xK, A^ = -27 (mK, \ y = 22 ^K. The 
RMS amplitude Qrms = 12.6 pK. Note that this value was 
calculated here without taking account of the existing corre- 
lation between the signal and noise in the COBE data, the 
inclusion of which would have given a sligh tly lower value 
(10.7 fJ,K), as explained in (Kogut et al. 1996). The principal 
vectors el, el, el point to the following positions in the sky 
(I, b): (335°, 2°), (250°, -63°), (64°, -27°). 

We have thus expressed the quadrupole signal as the 
sum of two orthogonal components proportional to the 
Legendre polinomial of second order, one of which is de- 
fined with respect to the minimum direction of the CMB 
quadrupole (±el) and the other with respect to its maxi- 
mum direction (±el). We can therefore take advantage of 
formulae (|^) and ([]) of Section 2. Although those equations 
were originally thought of in intensity terms, we can carry 



out our further analysis in terms of the brightness temper- 
ature by introducing temperature Stokes parameters, (T, 
Qt, Ut), bearing in mind that small changes in temperature 
parameters, e.g. AQt/T, are related to the corresponding 
changes in intensity parameters, i.e. AQ V /I V , via the fre- 
quency function 

d In I v xe x 



dlnT 



(10) 



where x = hvjkT is a dimensionless frequency. Since the 
change in the CMB in the direction of a galaxy cluster is very 
small, the resulting polarization degree will obviously have 
a frequency dependence defined by formula (p"o"|), P'(x) = 
f{x)(Q$ + U?) l ' 2 /T. 

After elementary calculations, which are essentially the 
conversion of the individual polarization signals that are pro- 
duced by each of the two orthogonal quadrupoles to common 
axes (see, e.g. Chandrasekhar 1950) and their subsequent 
summation, we obtain 



T'( 



T(9, 



r ! Ti + 0.9 



Ax I sin 2 9 cos 2 d> — — 



+A„ I sin 2 9 sin 2 c 



00 1 

)l+5> ,<n) 



y T (9,<j>) = 0.1r[Aa; (cos 2 #cos 2 (f> — sin 2 
+Aj,(cos 2 9 sin 2 <j> — cos 2 4> 



Ut(9,cI>) = 0.1t(X x - \ y ) cos 9 sin 20, 



P'(6,cf>,x) = 0.1 f^rTo 1 { [X x (l -sin 2 6>cos 2 



(12) 



(13) 



+ X y ( 1 — sin 2 9 sin 2 0)1 2 



4A I A H cos 2 6»j 1/2 , (14) 



where T(9, <f>) is the initial CMB temperature distribution, 
in which the main, isotropic component To = 2.73 K, and 
Tj and T n are the dipole and higher multipole harmonics. 
In the above, the Stokes parameters were defined so that 
if Qt > and Ut = 0, then the vector of polarization is 
orthogonal to 63. 

The polarization degree depends on the cluster po- 
sition on the celestial sphere in accordance with formula 
(p^). The angular dependence has two maxima, P max = 
0.1/(a;)rT _1 (X y — A^), in opposite directions, ±el. The 
maximal polarization signal expressed in absolute units is 
0.1(A y — X x )r = 4.9r fiK. In the directions of the min- 
ima and maxima of the CMB quadrupole, ±el and ±el, 
P = 0.1f(x)TT~ 1 X y and O.l/^rT^lAxI, respectively, and 
these four positions are saddle points in which the polar- 
ization is maximal in the plane 9 — n/2 and minimal 
in the perpendicular direction. There are also four posi- 
tions in the same 9 = 7r/2 plane for which P — 0: (f> — 
± arctan (— Xy/Xx) 1 ^ 2 and <fi = n ± arctan (— Xy/Xx) 1 ^ 2 . The 
geometry of the effect is sketched in Fig. 1. 

The sky-average (RMS) polarization induced by the 
CMB quadrupole can be found from expression ( |l4| ) to be 



X x Xy + Xy), 



(15) 
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U change their sign in these points. The area around a min- 
imum within which P/P ma x < 1/4 is a 14° radius circle for 
e 3 the adopted COBE values. 



maximum 




' Quadrupole \ 


minimum ' 




\p=o I 


| Polarization ^ 





The polarization caused by the CMB quadrupole is very 
small, reaching at the maximum P ma x ~ 2 ■ 10~ 6 /(a;)r or, 
equivalently, 0.1(r/0.02) /iK. This fact does not, however, 
automatically require the detector sensitivity threshold be 
accordingly low. Instead, one can derive the polarization cor- 
responding to a given direction by averaging polarization sig- 
nals taken from a large number of clusters located roughly 
around this direction (to be more precise, two opposite di- 
rections can be probed at the same time). In such a study, 
it is, in principle, even not necessary to have information 
on the optical depth of each particular cluster. Indeed, the 
polarization effect depends linearly on r, therefore it is the 
ensemble-average optical depth which is important. Further- 
more, all other possible polarization effects associated with 
clusters, which we discuss below, are not additive, in con- 
trast to the quadrupole effect, i.e. summing up signals from 
individual clusters will tend to smooth out their contribu- 
tion. 



Figure 1. The geometry of the polarization effect induced by the 
CMB quadrupole. The vectors ej, e"5 and ej define the eigensys- 
tem of the CMB quadrupole temperature anisotropy. The polar- 
ization effect has two broad maxima in the directions ±e$ orthog- 
onal to the plane which contains the minima (±el) and maxima 
(iej) of the quadrupole. In the same plane there are four direc- 
tions for which there is no polarization. 



or, using equation (]9J), 

Prms(x) = ^ f{x)T 9^1. (16) 

Note the y/6 enhancement factor in this expression. For the 
COBE data the sky-average polarization is 3.1r fiK, or 0.6 
of the maximal signal. 

Expanding expression ( |l4| ) in the vicinity of the maxi- 
mum directions, 9 — 0, 7T, gives 



AP 



\ x cos 2 (j) — Xy sin 2 (f> 



Xy X X 



(AO) 2 



(17) 



The quadratic dependence on the angular shift implies 
that the polarization maxima are very broad. Each of the 
two roughly elliptical areas around the maxima within 
which P j Pmax > 1/2 subtends about a quarter of the 
sky, as implied by the COBE measurement. Obviously, in 
these extended regions one has best chances to detect the 
quadrupole-induced polarizion. We note, however, that the 
measurement of the effect in the direction of its maxima is 
not sufficient for deducing all the 5 parameters that define 
the global quadrupole, so measurements in other directions 
are still necessary. 

In contrast to the maxima, the four directions in which 
there is no polarization signal are well resolved, as near them 



, y/—^xXy 
X v X x 



(18) 



4 OTHER POLARIZATION MECHANISMS IN 
CLUSTERS OF GALAXIES 

4.1 Polarization due to cluster transverse motion 



Sunyaev & Zel'dovich (1980) have shown that the proper 
motion of a galaxy cluster relative to the CMB will induce a 
polarization in the microwave signal measured from the clus- 
ter. They calculated this effect for the Rayleigh- Jeans region 
of the CMB spectrum and indicated (see formula (8) in their 
1980 paper) how to calculate the effect for any desired fre- 
quency. We now intend to find the frequency dependence of 
the produced polarization in explicit form. 

Let us assume that the initial microwave background 
is isotropic and has a black-body spectrum of temperature 
To. Now let an electron move through this background at a 
speed v. In the electron's rest frame, the radiation spectrum 
is also black-body, but with a temperature depending on 
the viewing direction. The spectral intensity in this frame is 
given by 



I u = C- 



(19) 



where C = 2(fcT ) 3 / \hc) 2 , = v/c, 7 = (1-/3 2 )" 1/2 , fj, is the 
cosine of the angle between the electron's velocity vector and 
the direction of incidence of a photon, and all the quantities 
are measured in the rest frame of the electron. The angular 
distribution of the background can be expanded in Legendre 
polynomials. Retaining only terms up to the second order 
in /? 2 , one gets 



I v = C 



e x - 1 



6(e* - l) 2 
e"(e" + l) o a2 ( 2 1 



(20) 



The minima are sharp because the Stokes parameters Q and 



Thus the quadrupole component of the initial distribution 
is 
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h 

Io 



e x {e x + 1) 
2(e* - l) 2 



2 . ,2 



(21) 



-* 0, corresponding to the Rayleigh- Jeans 
h/Io = B 2 . 

2 o2 



In the opposite limit, x — > oo, 



In the limit a; 
spectral region 

one has I2/I0 = l/2x 2 B 2 . The transformation of the result- 
ing polarization to the CMB frame will produce additional 
terms which are of higher than the second order in 3, so 
we will ignore them here. Using equations (|H|) and (^), we 
then finally obtain the polarization to be measured towards 
a cluster moving in the direction y! with respect to the ob- 
server, 



P' 



~2(e x 



(22) 



Here, Bt = 3(1 — /i ) ' is the transverse component of 
the peculiar velocity of the cluster. The polarization vector 
is perpendicular to the plane formed by the velocity vector 
and the observing direction. In the Rayleigh- Jeans frequency 
region, we readily find P' v = Q.1B 2 t, exactly the result of 
Sunyaev & Zel'dovich (1980). In the Wien part of the spec- 
trum (x — > 00), the degree of polarization of the measured 
signal is mu ch higher, a fact mentioned by Zel'dovich & Sun- 
yaev (1981). Fig. 2(a) shows the distribution of the polarized 
signal in the picture plane for a spherical cluster. 

It is easy to show that the frequency-integrated polar- 
ization J Q' v dv I J I' v dv — B 2 t, in agreement with the result 
of Audit & Simmons (EggS). 



4.2 Polarization effects due to finite optical depth 



Sunyaev & Zel'dovich (198C) have shown that there is an- 
other polarization mechanism associated with galaxy clus- 
ters which is connected with two consecutive scatterings of 
a photon within the intracluster gas, i.e. with finite optical 
depth of the gas cloud. Indeed, after first scattering in the 
intracluster medium, the CMB acquires an anisotropy due 
to the thermal and kinematic effects. The magnitude of this 
local anisotropic component is proportional to the product 
of r and r\ = kT e /m e c 2 for the thermal effect (T e being the 
temperature of the plasma) and 3 for the kinematic effect. 
Second scattering within the cluster will induce polarization 
signals of the order of r\T 2 and Bt 2 . In the first approxi- 
mation, the frequency dependence of the polarization will 
be the same as that of the thermal or kinematic effect. The 
resulting distribution of the degree and direction of polariza- 
tion over the cluster projected on the sky will be determined 
by the distribution of gas density and temperature within 
the cluster, on which depends the magnitude of the local 
quadrupole component. 

Before going further we would like to point out that 
the thermal and kinematic effects are not the only possible 
mechanisms capable of producing a local anisotropy. The 
local intensity distribution should also be distorted by grav- 
itational effects, the most important of which is the effect 
of a moving gravitational lense (connected with the cluster 
peculiar moti on). T he resulting polarization was calculated 
by Gibilisco ( 1997 ), and for very massive clusters it can be 
of comparable magnitude with t he effects we discuss. 

bunyaev & Zel'dovich ( |l98C| ) estimated the finite depth 
effects for a homogeneous spherical gas cloud. We will first 
check their results and then consider more realistic models 



\ • V \ 



(a) 



Y V. 



(b) 



(c) 



(d) 



(e) 



(0 



Vr-,/ 

>-f-\ 
' (g 



Figure 2. Predicted angular distributions of the polarized signal 
in the direction of galaxy clusters. In cases (a)-(d) the model gas 
cloud is spherical and an isothermal density law with n = 3/2 
is assumed, (a) The polarization oc /3 2 t caused by the peculiar 
motion in the indicated direction (see the insert in the center of 
the figure), (b) The finite optical depth polarization oc /3tr 2 for 
the same peculiar motion, (c) The finite optical depth polarization 
oc (kT e /m e c 2 )T 2 at frequencies x < 3.83. (d) The same but at 
x > 3.83. (e) The gas cloud is a spheroid with the ratio of the 
principal axes b/a = 0.8 and the symmetry (longer) axis in the 
X direction, (f) The same but the symmetry axis is inclined at 
45° to the picture plane, (g) The same but the symmetry axis 
points along the line of sight. The projected core of the cluster 
(as defined by the isothermal law) is shown by the dashed lines 
in all panels. 



for the gas distribution in clusters. It turned out that the 
simple formulae obtained in Section 2 are not readily appli- 
cable to the problem at hand. We have therefore chosen to 
follow a more straightforward approach, which is outlined 
below. 

We introduce a coordinate system in which the OZ axis 
is pointed from the cluster in our direction and the picture 
plane contains the OX and OY axes. Our task is to calcu- 
late the polarization that would be measured from a given 
(X, Y) projection point. For a given model we must specify a 
density law, p(X, Y, Z), where we, for simplicity, by density 
mean the number density of electrons. The gas temperature 
rj is assumed constant throughout the cluster in all of our 
models. 

The calculation runs as follows. To the effect from the 
(X , Y) point contribute all points (X, Y, Z) on the line of 
sight passing through this projection point. Therefore, one 
needs to find for each such point the angular intensity 
distribution (in which the isotropic component can be ig- 
nored) resulted from the first scattering, AI(X, Y, Z, 8, <f>). 
Here, the 8 and <f> angles define the viewing direction I: 
(sin 8 cos 4>, sin 8 sin (f>, cos 6). If we define the Stokes param- 
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eters with respect to the OX and OY axes (Q > 0, U — 
corresponds to a polarization in the OY direction), then the 
resulting polarization will be 



Ql(X,Y) = 



3<tt 



16tt 

x / dZp(X,Y,Z) / dCism 2 8cos(24>)AI„(X,Y,Z,8,<j>)[23) 



Ul(X,Y) 



3ctt 
16^ 



x / dZp(X,Y,Z) J dVism 6 sin {2(f)) AI V {X,Y,Z, 9,4)) ,(2A) 

where err is the Thomson scattering cross section. 

Equations (^) and ( pi] ) take account of the Rayleigh 
scattering function and of the fact that the polarization vec- 
tor of a photon scattered from the direction I towards the 
observer makes an angle (j> with the OY axis, so the rota- 
tion to the standard axes needs to be implemented (see, e.g. 
Chandrasekhar 1950). 

The intensity A/„ is equal to the product of the optical 
depth r as seen from the point (X, Y, Z) in the direction I 
and a geometry-independent factor which is determined by 
whether the kinematic or thermal effect has caused the local 
anisotropy. For the thermal effect 



AI v (X,Y,Z,0, 



■ T{X,Y,Z,0,4S)riIuh{x), 



where I v is the CMB spectral intensity, and 
xe x ( e x + 1 



Mx) 



1 V e x - 1 / 



(25) 



(26) 



(sunyaev & Zel'dovich 1972). In the case of the kinematic 
effect, 



m v {x, y, z, e, 0) = - T (x, y, z, e, d>)(i p)hf(x), 



(27) 



where f3 = /3(sin a cos V>, sin a cos ip, cos a) is the peculiar ve- 
locity vector, and the frequency dependence for the kine- 
matic effect / is the same as that for the polarization pro- 
duced by the CMB quadrupole, given by formula (|l^). Now 
we proceed to some examples. 



4-2.1 Homogeneous spherical cloud 
p(r) = 



po, 
0. 



r < 1; 
t > 1. 



For the f3r 2 effect the calculations in equations (p3|), 
( pi| ) can be completed analytically, yielding: 

Ql = ^(Xcos^-Fsin^v/l- (X* + f(x), (28) 

U' v = ^(XsinV+y cosV) V 1 - (X 2 + Y*)/3tT$I v f(x), (29) 

where we have introduced to = 2<ttPo, the central optical 
depth. 

The de rived result is identical to that of Sunyaev & 
Zel'dovich (1980). One important property of the effect de- 
is that its amplitude 



scribed by relations ( |28[ ) and (fc 

{(Q 2 +U 2 ) 1 / 2 ) depends only on the projected distance r from 
the cluster center. The azimuthal angle determines only the 
direction of polarization; the sign of the effect is different in 



the leading and trailing parts of the cluster (with respect to 
the direction of motion), as shown in Fig. 2(b). The maxi- 
mal polarization degree is reached at r max = l/v2) an d is 
equal to 0.025/3tTo f(x). In the Rayleigh- Jeans region, this 
is simply 0.025/3tTQ. One should note that in this case, in 
contrast to the effect of the order of f3fr, there is no net 
polarization from the cluster. 

In the case of the ryr 2 effect, one can easily verify that 
Ul(X, Y) = 0. The integration for Q' V (X, Y) in formula (||) 
can be implemented numerically. The polarization direction 
depends on the frequency of measurement in this case. At 
low frequencies, the polarization is radially directed in the 
picture plane - see Fig. 2(c). At x = 3.83, the thermal ef- 
fect changes its sign and at frequencies higher than this 
critical value, the polarization vector is orthogonal to the 
radius- vector, as shown in Fig. 2(d). The maximal polar- 
ization degree of 0.014?7Tq fr(x) is reached at r = 0.85. In 
the Rayleigh- Jeans spectral region, P = 0.028r?r 2 . A gain, 
like for the /3r 2 effect, the angular- integrated polarization is 
zero. 



4-2.2 Spherical cloud, King's density law 

p(r) = [l + (r/r c ) 2 }-". 

Fig. 3 shows how the polarization degree depends on 
the angular distance (r/r c ) from the cluster center for two 
values of n: 1 and 2. The maximal polarization is reduced 
by a factor of ~ 2 compared to the case of a homogeneous 
sphere for a given central optical depth To. 



4-2.3 Ellipsoidal cloud, King's density law 

It is known from observations that the shape of clusters of 
galaxies often significantly differs from spherical. The finite 
optical depth polarization effects will be modified notably in 
such situations. Let us consider an idealized model cluster 
which has a shape of the ellipsoid of revolution and whose 
density distribution is described by a modified isothermal 
law: p(r) = [1 + [X/af + (Y 2 + Z 2 )/b 2 ]- n . To make the 
situation more complete, we can rotate this ellipsoid around 
the OY axis. Examples of the resulting polarization pattern 
are presented in Fig. 2(e)-(g). The typical polarization de- 
gree is, in this case, of roughly the same magnitude as it 
was for the spherical cloud model, but its distribution in 
the picture plane can be absolutely different. In particular, 
generally, there is a polarized radiation in the direction of 
the cluster center - for example in Fig. 2(e) the degree of 
ellipsoidality (b/a) is only 0.8 but we see a maximum of 
polarization near the center, in sharp contrast to Fig. 2(c). 
Also important is that the net polarization from the cluster 
is no longer zero. 

It is clear from the above that the angular distribution 
of polarization caused by the finite depth effects for a given 
cluster will depend strongly on its shape and orientation, as 
well as on the distribution of gas within the cluster. 



5 DISCUSSION 

Microwave polarization measurements in directions of 
galaxy clusters provide a novel method for determining the 
quadrupole component in the CMB angular distribution, 
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Figure 3. Various polarization effects as functions of projected 
distance from the cluster center, calculated for the spherical 
isothermal density model (for n = 1 and 2). Shown is the fac- 
tor to be multiplied by: (3 2 to (solid lines; for this effect as well 
for the CMB quadrupole effect the angular dependence simply 
follows the angular dependence of the line of sight optical depth), 
/3t7~Q (dashed lines), and (kT e /m e c 2 )Tg (dash-dotted lines), where 
to is the central optical depth. 



which does not require mapping of the whole sky, as is the 
case with the conventional method of measuring the CMB 
anisotropies. What is more important is that this method 
allows one to measure the CMB quadrupole component as 
would be seen by observers that are placed on distant clus- 
ters, including those at z ~ 0.5-3. This information, which is 
otherwise unaccessible, has significant implications for cos- 
mology, which should raise interest of experimentators to 
microwave polarization measurements of distant clusters. 
Such measurements can at the same time provide another 
piece of information valuable for cosmology, namely that 



on the evolution of the p eculiar velocities of clusters (3un 



yaev & Zel'dovich 1981). The most serious concern is, of 



course, the small magnitude of the polarization induced by 
the CMB quadrupole, which is expected to be less than 
~ 0.1(t/0.02) fiK. However, as we mentioned in Section 3, 
one can overcome the detector sensitivity limit by probing 
a large number of clusters. 

Another problem is the presence of several other po- 
larization effects connected with clusters of galaxies. These 
effects, which are interesting for themselves because they 
could, in principle, provide us such important information 
as on cluster tangential velocities and the distribution of in- 
tracluster gas, can produce a polarization comparable in size 
with that due to the CMB quadrupole, and will thus be a 
source of significant noise for observations aiming at the de- 
tection of the quadrupole-induced signal. Fortunately, all of 
these effects have a tendency to vanish in the course of aver- 



Figure 4. Various polarization effects as functions of frequency: 
the CMB quadrupole effect and the /3|t effect (solid line), the 
/3tr 2 effect (dashed line), and the (kT e / m e c 2 )r 2 effect (dash- 
dotted line). The latest effect changes its sign at x = 3.83. 



aging over a large cluster sample. Let us see how significantly 
the different polarization mechanisms could contribute to an 
individual measurement. 

The polarizations due to the three basic effects are of 
the order of 0.1/3 t 2 r, 0.01/3 t T 2 and O.OI77T 2 in the Rayleigh- 
Jeans spectral region (the coefficients at the last two effects 
are model-dependent, but typically not differ much from 
the quoted values). The first effect becomes equal to the 
CMB quadrupole polarization effect at its maximum for the 
tangential velocity ~ 1300 km/s (as implied by the COBE 
data). This value is considerably higher than the average 
peculiar velocity expected for clusters of galaxies (see, e.g. 
Bahcall & Oh 1996). Moreover, the sky-average of the po- 
larization signal due to the CMB quadrupole is ~ 60% of 
the maximum value, whereas the velocity effect is propor- 
tional to the square of /3 t , and therefore for typical values, 
Pt ~ 400 km/s, the polarization produced by the CMB 
quadrupole should be several times stronger in the Rayleigh- 
Jeans region of the CMB spectrum. Among the considered 
effects, only the O.Olr^r 2 polarization might compete with 
the quadrupole-induced polarization for rich clusters. 

The above values were calculated for the Ray leigh- Jeans 
spectral region, which corresponds to centimeter and mil- 
limeter wavelengths. However, the polarization degree will 
strongly depend on frequency if one performs measurements 
in the short millimeter and, especially, in the submillime- 
ter band. The frequency dependencies for the different ef- 
fects are shown in Fig. 4. Note that the polarization due to 
the CMB quadrupole has the same frequency dependence 
as the f3 t r 2 polarization, both growing linearly (~ x) in the 
Wien part of the spectrum. The other two effects rise more 
steeply, ~ x 2 , in the Wien region. The different frequency 
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behaviours provide the opportunity to separate the signals 
of different nature by carrying out multi-band observations. 
On the other hand, it is clear that there is more hope to 
detect the polarization due to the CMB quadrupole in the 
Rayleigh- Jeans region, where its contribution is relatively 
more significant. 

The effects being discussed also depend in different ways 
on the projected angular distance from the cluster center. 
The CMB quadrupole effect and the /3fr effect are propor- 
tional to the optical depth through a given line of sight, 
and therefore both lead to a finite polarization of the sig- 
nal integrated over the picture plane. The other two effects 
have a tendency to vanish in the net signal. This fact sig- 
nificantly enhances our chances to determine both the CMB 
quadrupole and the tangential velocities of clusters by means 
of polarization measurements - one just has to accumulate 
the microwave signal from the whole cluster (or at least 
from an area which is symmetric around its center). Fur- 
thermore, the 0tT polarization will vanish upon averaging 
over a large set of clusters, unless the directions of their 
tangential motions are correlated, enabling the extraction 
of the pure quadrupole-induced signal. Fig. 5 demonstrates 
the complexity of the overall problem. It shows how the po- 
larization pattern shaped purely by the CMB quadrupole 
anisotropy will be modified by the other polarization effects 
and undergo further changes if one tunes to a different mea- 
surement frequency. 

ft is worth noting that other polarization effects, apart 
from those discussed in this paper, may appear strong under 
some circumstances. For example, the signal from a clus- 
ter can be polarized if there is a strong compa ct source of 
unpola rized rad i o emission insi de the cluster (Sunyaev & 
Zel'dov fch 198C| ; Sunyaev 1982). The polarization pattern 
in this case will be similar to that produced by the r/r 2 ef- 
fect - see Fig. 2(d). Of course, the frequency dependence 
in this case, which should be power-law like, will markedly 
distinguish itself. Another way to gi ve rise to pola rization 
is through gravitational mechanisms ( Gibilisco 1997 ), which 
we alre ady mentioned above. Finally, there is a polariza- 



tion wl ich is intrinsic to the CMB. This polarization sig- 
nal, which together with the intensity fluctuations has been 
brought to us largely from the epoch of decoupling but may 
also bear imprints from later epochs, is expected to be very 
small at sub-arcminute scales which correspond to the typi- 
cal dimensions of cluster cores (see, e.g. Seljak & Zaldarriaga 
1998; White 1998). ft should also be noted that the typical 
Faraday rotation measures of galaxy clusters are too small 
to produce a significant depolarization of the discussed ef- 
fects for wavel engths below ~ 10 cm, i.e. in the doma in of 
interest to us (Lawler & Dennison 1982; [garazin 1988 ). 
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Figure 5. Combined effect of the CMB quadrupole and other 
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in the direction of the maximum of the quadrupole-induced po- 
larization, has a transverse velocity of 600 km/s (/3t = 0.02), 
central optical depth r = 0.01 and gas temperature 5 keV 
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size of the circle shown on each panel, (a) Only the quadrupole- 
induced polarization, (b) The /3 2 r effect has been added, (c) The 
effect (kT e /m e c 2 )T 2 has been added. Cases (a)-(c) are for the 
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